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I.  INTRODUCTION 

The  problem  of  electromagnetic  radiation  and  scattering  from 
perfectly  conducting  bodies  of  revolution  of  arbitrary  shape  has  been 
studied  by  Mautz  and  Harrington  [1].  An  integro-dif ferential  equation 
was  obtained  from  the  potential  integrals  plus  boundary  conditions  at 
the  body  and  was  solved  by  the  method  of  moments.  The  simplification 
derived  from  a  rotationally  symmetric  body  is  that  all  physical 
quantities  can  be  expanded  in  Fourier  series  in  the  azimuth  coordinate 
and  quantities  of  different  harmonics  can  be  treated  independently.  A 
slightly  different  formulation  based  on  the  reaction  theorem  was  also 
used  by  Chuang  [2]  to  solve  the  same  type  of  radiation  and  scattering 
problem. 

In  this  report,  the  same  formulation  given  in  [2]  is  employed  to 
set  up  generalized  network  equations  for  perfectly  conducting  bodies  of 
revolution  for  application  to  the  radiation  problem  of  conical-type  horn 
antennas.  These  antennas  include  smooth-wall  conical  horns,  corrugated 
horns  and  dual-mode  horns.  The  corrugated  horn  pattern  is  simulated  by 
the  H-plane  pattern  of  a  conical  horn  with  the  same  inner  dimension  as 
the  corrugated  horn  and  is  assumed  to  be  circularly  symmetric.  Several 
examples  are  shown  and  compared  with  measured  patterns  in  order  to 
validate  this  approach.  Appendix  A  shows  the  detailed  derivation  of  the 
elements  associated  with  the  impedance  matrix.  The  Body  of  Revolution 
Code  has  been  modified  and  implemented  into  the  OSU  Reflector  Antenna 
Code.  One  can  look  at  the  BR:  Command  of  the  User's  Manual  of  the 
Reflector  Antenna  Code  for  more  detail. 
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II.  REACTION  FORMULATION 


Consider  a  primary  source  of  an  electric  current  J.  illuminating  a 
perfectly  conducting  horn  type  scatterer.  The  scatterer  occupies  a 
volume  T  enclosed  by  a  surface  S  as  shown  in  Figure  2.1(a).  Currents 
are  induced  on  the  scatterer  which  in  turn  will  radiate  such  that  the 
total  fields  are 


E  =  , 

(1) 

H  =  +  H®  , 

—i  —i  — s  —s 

where  E  ,  H  and  E  ,  H  are  the  primary  and  scattered  fields, 
respectively. 

Now  let  us  assume  that  the  horn  walls  have  zero  thickness  and 
replace  the  scatterer  with  free  space  as  shown  in  Figure  2, 1(b).  From 
the  equivalence  principle,  one  can  use  the  equivalent  surface  current  S 

s 

to  find  the  total  scattered  field.  The  surface  current  ^  on  the 

s 

surface  S  is  given  by 


=  n  X  (iTj  - 


(2) 


while  the  magnetic  surface  current  vanishes  since  the  scatterer  is  a 
perfect  conductor.  In  Equation  (2),  n  is  the  unit  outward  normal  to  the 
surface  and  and  ^2  the  total  scattered  fields  outside  and  inside 
the  horn.  Further  let  us  introduce  an  arbitrary  current  into  region 


2 


(a)  original  problem 


A 


Jg  =  n  X  (  H ,  -  H  2  ) 


(b)  free  space  equivalence 
Figure  2.1.  Geometry  of  horn  scattering  problem. 


T  originally  occupied  by  the  scatterer.  Assume  that  ?  will  radiate 

m 

fields  if  into  free  space.  Then  from  the  reciprocity  theorem  we 
arrive  at  the  reaction  equation 


Is  •  It 

■  -  It 

=  -  f  E*"  •  1  d\) 
Jv. 


(3) 


Note  that  the  volume  integral  is  performed  over  the  region  T  as  shown  in 
Figure  2.1(a)  although  the  thickness  of  the  horn  surface  is  assumed 
zero.  The  condition  ?  =  0  on  the  surface  has  been  used  in 

Equation  (3).  The  last  integral  in  Equation  (3)  is  performed  over  the 
region  where  the  source  is  located. 

Let  the  z-axis  be  the  axis  of  revolution  and  t  the  distance 
measured  along  a  profile  line  of  the  conducting  body  as  shown  in  Figure 
2.2.  Because  of  the  orthogonality  of  Fourier  series,  different 
harmonics  in  azimuth  coordinate  ((J>)  can  be  considered  separately.  For 
each  harmonic,  the  induced  currents  on  the  conducting  body  can  be 
written  as 


J 

s 


^cos  m<fi' 

^sin  m<()  ' 

=  t  f(t)< 

Isin  m<f>; 

-  +  <l>g(t)< 

,-cos  m(J)> 

(4) 
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Figure  2.2.  Coordinate  system  for  horn  scattering  calculation  by  moment 
method. 

where  m  can  be  zero  or  any  positive  integer.  For  a  conical  type  horn 
antenna,  we  choose  a  short  electric  current  dipole  oriented  in  the  x- 
direction  and  located  on  the  z-axis  inside  the  horn  as  the  primary 
source.  Therefore  only  the  upper  sinusoidal  function  in  (4)  will  be 
used  with  m=l. 

Following  standard  procedures  of  moment  methods,  Equation  (4)  is 
then  expanded  in  a  finite  series 


-  Q  (t)  .  . 

t  I  v  /— V-  cos  m<|>  +  <(»IT  (t)  sin  m* 
n  Rp(t;  n  n 


(5) 


where  p(t)  is  the  distance  from  a  point  t  on  the  surface  to  the  axis  of 
revolution  as  shown  in  Figure  2.2.  Qj^(t)  and  Pjj(t)  are  expansion 
functions  defined  as 


Qn(t) 


sin  k(t-t  1 ) 
_  n-l*^ 

sin  k(t  -t  . 
'  n  n-1 


)  ’ 


if  t  >t>t 
n-  ~ 


n-1 


sin  k(t_^^j-t) 


if  t  <t<t  . 
n-  -  n+1 


=  0 


otherwise, 


(6) 


and 


P 

n 


(t) 


1  , 


=  0  , 


if 


t  <t<t  . 
n-  -  n+1 


otherwise, 


(7) 


where  t^  are  sampling  points  along  the  generating  profile  line  of  the 
body  of  revolution  as  shown  in  Figure  2.3(a). 

If  the  test  sources  are  chosen  to  be  the  same  as  the  expansion 
currents.  Equation  (3)  then  reduces  to  a  matrix  equation. 


[z"]  [z'*] 

Tv']- 

Iv^]. 
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The  elements  of  the  various  matrices  are 


-tl 

t  cos  m<}>  5 — 7—^  •  E  ds 
kp(t) 


sin  m<|)  *  E  ds 


Q_(t) 

t  cos  .  E  ds 


<j)  sin  m<|)  ds 


t  r  -  -i 

Vi  =  t  cos  mi>  •  E  ds 


=  ij)  sin  m<f>  Pj^(t)  •  ds 


with  m=l, 
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E  =  -E'  ON  S 


Figure  2 


A 


(a)  scattered  field  produced  by  horn  wall  currents 


.3. 


(b)  field  from  single  expansion  current  ij  t 
Expansion  of  the  surface  current  ?  • 


Note  that  and  are  electric  fields  generated  by  electric  current 
modes  and  1^,  respectively.  From  the  reciprocity  theorem  it  can  be 
shown  that  the  impedance  matrix  is  symmetric. 


A.  The  Impedance  Matrix 

In  this  section  the  impedance  elements  will  be  discussed  in  detail. 
From  Equation  (9)  the  impedance  elements  have  the  general  form: 


ds 


=  j» 

77 

ds 

=  j(0 

4 

(v 

yi 

1^7  •  A^jj  ds 

. ^  -jkR 

(10) 


The  derivation  of  this  equation  is  given  in  Appendix  A  (in  detail).  In 
view  of  the  following  formulas, 


7  •  J 


1 

P 


_9 

3t 


(11) 
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t*t'  =  sin^  sin5'  cos(<(>-<|)' )  +  cos^  cosC'  > 

=  sin^  sin(<j)-<|)' )  , 

(}).  (ji'  =  cos(<|>-4>')  (12) 

where  ^  is  the  angle  between  t  and  z,  the  various  impedance  elements 
become 


,tt 

■'In  "  ,,  2 

4k 


sin^  sin^' (G^^j^+G^_j)+2  cos^  cos^'G^  Qj^(t)Q^(t') 


•d  t  •  d  t '  , 


4k2 


{kpsin^'(G^^j-G^_l)P^(t)Qi(t')-^G^Pn(t)Qi(to}dt.dtS 


4k2 


{kp'sin^(Gm^,-G^_l)Qn(t)Pi(t')-^  )]dt-dV  , 


AA  J 

_ O 

^In  ■  2 

4k 


7<(>t  „t<l> 

^In  =  nl  ' 
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(13) 


These  equations  are  derived  in  detail  in  Appendix  A. 

The  variables  in  Equations  (13)  are  defined  as  follows: 


G  = 
m 


-jkR 


cos  m4> 


d<fr 


(14) 


R  = 


(p-p)  +  4pp' 


• 


(15) 


The  evaluation  of  has  been  discussed  in  detail  by  Mautz  and 
Harrington  [1]. 


B.  Excitation  Column 

The  primary  source  is  chosen  to  be  an  x-directed  on-axis 
electric  dipole  with  dipole  strength  The  excitation  element  vj 

associated  with  one  of  the  monopoles  of  the  current  t  is  given  by 


ds' 


in  which  is  the  electric  field  radiated  by  the  primary  source  1  and 
the  integration  is  carried  over  a  ring  of  surface  current.  From 
reciprocity,  one  can  calculate  by 
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where  is  the  electric  field  radiated  by  the  current  mode  t  and  the 
integration  is  carried  over  the  space  at  which  the  primary  source  is 
located.  Since  the  primary  source  is  an  on-axis  dipole,  one  can  express 
as 


=  X  P^8(x)8(y)5(z-z^) 

where  8(x)8(y)8(z-z^)  indicates  the  dipole  is  located  at  (0,0, z^). 
Thus,  is  calculated  by  the  following  equation 

Vj  =  X  Px  •  f^(0,0,z^)  . 

The  geometry  for  the  calculation  of  f^(0,0,z^)  is  shown  in  Figure  2.4. 
The  monopole  is  located  between  t=tj  and  t=t2  with  P  being  the  field 
point.  The  current  distribution  on  the  monopole  is  given  by 


t  j:  =  t 


IjSin 


kp(t) 


jk(t2-t)  +  I2  sin  k(t-tj^)j 


sin  kd 


cos 


since  piece-wise  sinusoidal  functions  have  been  used  to  expand  the 
surface  current.  The  electric  field  fj^(0,0,z^)  radiated  by  t  jJ  can  be 
found  by  the  following  equation 


where  is  the  magnetic  vector  potential  and  §  is  the  scalar  potential. 
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Figure  2.4.  Geometry  for  the  radiation  of  a  monopole  with  end  points 
t^ ,  t„  and  the  field  point  0  on  the  axis  of  revolution. 


The  vector  potential  is  given  by 


^  f  *  t  e-j*^ 

^1= 


^  Ij^sin  |^k.(t2-t' )]  +  I2  sin  l^k(t'-tj^)  j 


0  Jt, 


kp(t') 


sin  kd 


cos  m4> 


p(t')  dt'  d(|> 


-  p2jl  frt,  .  I-sin  k(t„-t')  +  sin  k(t'-t.)  -jkR  " 

A,  =  t-  cosm<f«  t  - : — - .  _p  dt'  >d<t> 

1  k  ^  ^  sin  kd  4itR 

JU 


1  -> 
r  cosmt})  A'  d<() 
R  0  J- 


From  the  rotational  symmetry  of  the  problem,  one  can  calculate  the 
electric  field  associated  with  and  then  integrate  over  i>  to  find 
The  electric  field  radiated  by 


.1-  sin  k(t^-t')  +  sin  k(t'-t-) 

?(t)  .  t  Kt)  -  t  - 
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is  similar  to  the  radiation  of  a  z-directed  current  element  with  the 
gamp  current  distribution.  The  cylindrical  components  of  the  field  are 
given  as  [3] 


E 


t  4ltj  sin  kd 


-i“2 

-jkR-i' 

f  \  G  ^ 

[[ij  -  cos  kdj  g,  + 

(ij  -  cos  kd]  J 

and 


E  = 
n 


4nj 


p  sin  kd 


(^1 


-jkR. 


-jkR-} 


j  sin  kd 


-jkR. 

+  (Ijcos  kd  -  l2)e  cos  0^^  +  (I2C0S  kd  -  Ij^)e  cos  62 

where  h  is  the  free  space  wave  impedance  and  E^  is  the  component  along 
the  axis  of  the  monopole  and  E^  is  the  component  orthogonal  to  the  axis 
The  electric  field  is  thus  given  by 


= 


t  E^  +  n  E 


z  >  z 
o  c 


or 


*1  -  t  Ef  .  n’  E„ 


t  E^  -  n  E 
t  n 


z  <  z 
o  c 


as  shown  in  Figure  2.5. 
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From  Figure  2.2, 


t  X  p  =  sin^  +  z  cos^j  X  [-PP^^  +  2(z^-Zj)j 
=  |^-(z^-Zj^)sin^  -  pj^cos^j<j> 

=  -  [^(z^-Zj^)sin^  +  p^cos^ji|)  . 

When  z^  >  z^,  t  x  p  is  in  the  -  <|)  direction  so  that  |^(z^-Zj^)sin^+Pj^cost 
>  0;  when  z^  <  z^,  t  x  p  is  in  the  <|)  direction  so  that  ^(z^-z^)sin^  + 
pj^costj  <  0.  Consequently,  the  electric  field  can  be  written  as 

V 

E'  =  t  +  sgn[(z^-Zj)sin^  +  p^cos^]  n  E^ 

or 

=  (p  sin^  +  z  cos^E^  +  sgn[ (z^-Zj^)sin^  +  pj^cos^]  • 

(-P  cos^  +  z  sin^)  E^ 


or 
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%  =  p|^sin^  -  sgn[  (z^-Zj^)sin^  +  p^cos^]  cos^ 

+  z<^cos^  E^  +  sgn[ (z^-z^)sin^  +  pj^cos^]  sin^  E^j> 

=  |x  cos<^  +  y  sinifijl^sin^  E^  -  sgn[ (z^-Zj^)sin^  +  pj^cos^]  cos^  E^^ 
+  z|cos^  E^  +  sgn[ (z^-Zj^)sin^  +  p^cos^]  sin^  E^| 


Consequently, 


riK 


?,(0,0,z  )  =  ^  cos  m<()  ^'(0,0, z)  d(() 
1  O  K.  ^  1 


or 


fj^(0,0,z^)  = 


0  ,  m>l 

X  ^^sin^  E^  +  sgn[  (z^-Zj^)sinJ.  +  pj,cos^]  cos^  E^j>,  m=l 
z  ^^cos^  E^  +  sgn[  (z^-Zj^)sin^  +  p^cos^Jsin^  E^j',  m=0 


t  . 


Thus,  the  excitation  element  is  given  by 


vi  =  X  P  •  f-,(0,0,z  ) 
1  X  1^  ’  ’  o 


=  Px  E  {®t  “  sgn[(z^-Zj)sin^  +  p^cos5]E^  cos^j'  ,  m=l 


(17) 
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In  Equation  (17),  is  the  dipole  strength  of  the  primary  source,  the 
function  sgn(x)  is  defined  as 

sgn(x)  =  ±  1,  if  X  ^  0, 

and  (0,z  )  and  are  the  coordinates  of  the  point  P  and  t^^. 

The  excitation  element  associated  with  the  current  mode  <(>  is 
obtained  accordingly  by  the  equation 


<|.  jJ  •  ds 


=  -J» 


1 


(v  •  jJ) [v  •  l^jjds 


(18) 


The  vector  potential  due  to  a  perpendicular  on-axis  dipole  is 


-jkR 

ye"’ 

:?i  r.  0 

A  =  X  P  — — 
X  4nR 


(19) 


A 

where  R  is  the  distance  from  the  on-axis  dipole  to  the  current  J£.  The 
current  is  given  by 

i  J*  =  i  [ij  Pj^(t)  +  P^^j(t)]sin  m«|.  . 


In  view  of  the  form  of  Jp  the  excitation  element  associated  with  the 
current  mode  Pj^(t)  sin  m<|>  is  discussed  next. 
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From  Equation  (18),  the  excitation  element  associated  with 


<f>  Pj^(t)  sin  m(j>  is  given  by 


JW 


{•f- 


sin  m(p 


i  _  L[^ 
L 


(<(• 


P^Ct) 


sinm 


y  • 


where 


y  e  ye 

•  P  =  4»-x  P^  -  =  -  sin<|)  P^  ° 


x  AnR 


■  X  4llR 


V  •  (<|)  Pj^(t)  sinm<(>)  =  ^  (P]L^ p  ’ 


and 


7  • 


X  P 

X 


AnR 


o 

X  An 


y  P 

o  X 
All 


aR 


aR 

ax 


_  [_  jk  -jkR 

'Art  I  R 


x-x' 

R 


For  the  on-axis  dipole,  x'=0  and  x=p  cos4>  so  that  7'A  becomes 
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.  _  llo!^  L  ik  _  1_ 

"  "  4rt  I  R 


p  cos<^  . 


Substituting  the  above  relations  into  V^,  the  following  equation  is 


obtained: 


sin<(i  sinin<|» 


jwu  pit  ro  nkl"^  '1 

Pj^(t)|k  sin<|)sinm<()  -  m  +  ^Jcosm<|>cosij>J 

» 0  *  t 


■pdtd<t> 


The  integration  over  ^  results  in 


2ii  p  ,  m^l 

sin<|>  sinni<|>  d+  =  - 
0  U  ,  m=l 


2ii  p  ,  m^l 

cos<(>  cosm<|>  d^  =  ' 

0  U  ,  ni=l 


For  the  horn  radiation  problem,  m=l,  and  becomes 


♦  _  ^::r2  p  p  (t)  fk^  _  lis  _  u 
1  -  4k^  X  .  *^1^^^  r  R  R'j 
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The  complete  is  then  the  sum  of  two  such  equations  (Equation  (20)) 
associated  with  the  two  current  modes  comprising  the  current 


C.  Far  Zone  Fields 

The  vector  potential  in  the  far-field  zone  due  to  a  current  mode  of 
Equation  (5)  is  given  by 


U  e 

i  =  ^ 


-jkr 


n  4Jir 


J  p'dt'dr 


-jkr 


4  nr 


t'  COS  I  — I — +  (})'  sin  m<{)  I  P  (t') 
n  kp'  ^  n  ^ 


^jk[p'sinecos((|)'-<)))+z'cose] 


(2.1) 


In  view  of  the  relations  (excluding  the  r  component) 


<}>'  =  ij)  cos(if)'-<j))  -  e  sin( ij)' -<j))  cos6  (22) 


t'  =  <|)sin^'  sin( <f>' -i}))  +  e[sin^'  cos6  cos(il>'-4>)  -  cos^'  sin9]  . 

(23) 


Equation  (21)  then  reduces  to 


A 

,<P  _  o 

nj> 


-jkr 


75 -  I  sin  m(j)  j 

2  r  n 


.m-1 


J'(kp'sin9)P  (t')e^^^  *^°^^p'dt' 
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-jkr 


ne 


2  r 


I  cos 
n 


m«|»  j 


.m-1 


Q  J  (kp'sin0)P 
kp'sinG  m'  '  ' 


n 


-jkr  T-t 

-u  e  I  -1 

.t  n  .  .  ,  .m-1 

'"n*-  2  r  r  "*  J 


sin9)0_(  t '  t ' 


-jkr  _t 

.t  “o' 

*n9-  -rr-  r  ““  "♦  ^ 


.m-1 


[J^  (kp'sinG)  sin^'cos0 


-  j  cos^'  sin0  J^(kp'sin0)]Qj^(t' 


(24) 


where  J  is  a  Bessel  function  and  J'  is  its  derivative  with  respect  to 
m  m 

the  argument. 

For  the  radiated  field  of  a  conical  horn  antenna,  one  also  needs  to 
include  the  radiation  due  to  the  primary  dipole  source.  In  the  far- 
field  zone,  the  vector  potential  of  the  perpendicular  on-axis  dipole  of 
strength  is 


^  o  -  PqPjj  -jkr+jkz^cos0 

A  =  (0  COS0  COS'!)  -  <)>  sin<|>)  e 


(25) 


The  total  radiated  field  can  then  be  obtained  from  (24)  and  (25)  as 


r  «s  ■  -  d«'„  *e,  total 

i 

k  E.  =  -jdoy  A.  ^  ^  , 
o  <j),  total 
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III.  EXAMPLES 


In  this  section,  six  horn  antennas  are  used  as  examples  to  validate 
the  moment  method  computer  code.  The  calculated  patterns  are  compared 
with  the  measured  patterns. 

A.  Prime  Focus  Conical  Horn 

The  geometry  of  this  horn  antenna  is  shown  in  Figure  3.1.  The 
radiation  patterns  of  this  horn,  calculated  at  11  GHz,  are  shown  in 
Figure  3.2.  The  measured  patterns  are  given  in  Figure  3.3  for 
comparison. 

B.  Prime  Focus  Corrugated  Horn 

The  geometry  of  this  corrugated  horn  is  given  in  Figure  3.4.  The 
inner  dimension  of  this  horn  is  used  as  a  conical  horn  and  the  H-plane 
pattern  of  this  simulated  conical  horn  is  calculated  at  11.0  GHz.  The 
resulting  pattern,  which  is  given  in  Figure  3.5,  is  treated  as  the 
pattern  of  the  corrugated  horn  and  assumed  circularly  symmetric.  The 
measured  principal  plane  patterns  are  shown  in  Figure  3.6. 

C.  Prime  Focus  Dual -Mode  Horn 

The  geometry  of  this  horn  is  shown  in  Figure  3.7.  The  moment 
method  calculated  patterns  are  shown  in  Figure  3.8  with  the  measured 
patterns  shown  in  Figure  3.9. 
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D.  Conical  Horn  for  Cassegrain  Reflector 

The  geometry  for  this  horn  antenna  is  shown  in  Figure  3.10.  The 
radiation  patterns  of  this  horn,  calculated  at  11.0  GHz,  are  shown  in 
Figure  3.11.  The  measured  patterns  are  given  in  Figure  3.12  for 
comparison. 

E.  Corrugated  Horn  for  Cassegrain  Reflector 

The  geometry  of  this  corrugated  horn  is  given  in  Figure  3.13.  The 
inner  dimension  of  this  horn  is  used  as  a  conical  horn  and  the  H-plane 
pattern  of  this  simulated  conical  horn  is  calculated  at  11.0  GHz.  The 
resulting  pattern,  which  is  given  in  Figure  3.14,  is  treated  as  the 
pattern  of  the  corrugated  horn  and  assumed  circularly  symmetric.  The 
measured  principal  plane  patterns  are  shown  in  Figure  3.15. 

P.  Dual-Mode  Horn  for  Cassegrain  Reflector 

The  geometry  of  this  horn  is  shown  in  Figure  3.16.  The  moment 
method  calculated  patterns  are  shown  in  Figure  3.17  with  the  measured 
patterns  shown  in  Figure  3.18. 
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Con  i  c  a  1  l-l  o  r  n  for  P  r  i  me  Focus  fi  n  L  e  n  n  a 


Figure  3* 1(a).  Geometry  of  the  prime  focus  conical  horn. 


Figure  3.1(b).  Equivalent  geometry  of  the  prime  focus  conical  horn  for 
moment  method  calculation. 


FF  PRTTERN  (DB)  .  FF  PATTERN  (DB) 


Cfll.CULflTEO  PRIME  FOCUS  CONICAL  HORN  CO-POL 
PHI  (DEC)  -  0.00  FRE(3  (GHZ)  ~  11.00 


.  -ISO.  -120.  -90.  -60.  -30.  0.  30.  60.  90.  120.  150.  180. 


THETFI  (DEGREES) 
(a)  H-plane 


(b)  E-plane 

Figure  3.2,  Calculated  patterns  of  the  prime  focus  conical  horn  at  11.0 
GHz. 
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Figure  3.4(b).  Equivalent  conical  horn  of  the  prime  focus  corrugated 
horn  for  moment  method  calculation. 


FF  PATTERN  (08) 


3.5. 


Calculated  pati 
11.0  GHz  (simul 


PFCR_OCQ_l  lOO.OflT  PF  CORRUCflTEO  FEED  CO-POL  -180/180  PHI  =  0. 
PHI  lOEC)  =  0.00  FREQ  (GHZ)  =  11.00 


PFCR_90C0_1 lOO.DRT  PF  CORRUGRTED  FEEO  CO-POL  -180/180  PHI=90. 

PHI  lOEG)  =  90.00  FREQ  (GHZ)  =  11.00 


(b)  E-plane 


Figure  3.6.  Measured  patterns  of  the  prime  focus  corrugated  horn  at 
11.0  GHz. 
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Dual  Mode  I  lorn  for  Prime  Focus  DnLenna 


Phase  Center 


/  \  /  la 

/  \  “  S 

J-'-l  :  I 


- 0.519“ 


Figure  3.7(a).  Geometry  of  the  prime  focus  dual-mode  horn. 
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Figure  3.7(b).  Equivalent  geometry  of  the  prime  focus  dual-mode  horn 
for  moment  method  calculation. 
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PF0M_0C0_1100.DflT  PF  OUflL-MODE  FEED  CO-POL  -180/180 
PHI  (DEG)  =  0.00  FREQ  (GHZ)  =  11,00 
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-90.  -60.  -30,  0.  30.  60.  90, 

THETfi  (DEGREES) 

(a)  H-plane 


PFDM_90C0_i 100.  CRT  PF  OUflL-MODE  FEED  CO-POL  -180/180 
PHI  (OEG)  =  90,00  FREQ  (GHZ)  =  11.00 


-1^.  -150.  -120.  -90.  -60.  -30.  0.  30.  60.  90.  120. 

THETfl  (DEGREES) 

(b)  E-plane 

Figure  3.9.  Measured  patterns  of  the  prime  focus  dual 
GHz. 
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Figure  3.12.  Measured  patterns  of  the  conical  horn  for  Cassegrg 
reflector  at  11.0  GHz. 
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Corrugated  H orn  For  Cassegrain  Cintenna 


Figure  3.13(a).  Geometry  of  the  corrugated  horn  for  Cassegrain 
reflector. 


Figure  3.13(b).  Equivalent  conical  horn  of  the  corrugated  horn  for 
Cassegrain  reflector. 
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Figure  3.14.  Calculated  patterns  of  the  corrugated  horn  for  Cassegra 
reflector  at  11.0  GHz. 
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Dual  Mode  Horn  for  Cassegrain  Antenna 


Figure  3.16(a).  Geometry  of  the  dual-mode  horn  for  Cassegrain 
reflector. 


Figure  3.16(b).  Equivalent  geometry  of  the  dual-mode  horn  for 
Cassegrain  reflector. 
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FF  PATTERN  (DB)  FF  PATTERN  (DB) 


FF  PATTERN  (OB)  .  FF  PATTERN  (DB) 


IV.  SUMMARY 


The  method  of  moments  is  used  in  this  report  to  calculate  the 
radiation  patterns  of  horn  antennas  which  can  be  treated  as  bodies  of 
revolution.  The  simplification  derived  from  a  rotationally  symmetric 
body  is  that  all  physical  quantities  are  expanded  in  Fourier  series  in 
the  azimuth  coordinate  and  quantities  of  different  harmonics  are  treated 
independently  [1].  The  harmonics  related  to  the  horn  antenna  problem 
are  those  with  m=l,  i.e.,  cos4>  and  sin<j>. 

Several  examples  are  shown  and  compared  with  the  measured  patterns. 
Very  good  agreement  has  been  obtained  between  calculations  and 
measurements.  One  of  the  advantages  of  using  the  moment  method  instead 
of  the  well-known  aperture  Integration  method  for  the  horn  antenna 
problem  is  that  the  backlobes  of  the  radiation  patterns  can  be 
calculated.  Also  the  aperture  integration  is  not  accurate  for  very 
small  aperture  antennas  with  aperture  diameters  less  than  1.2X. 
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APPENDIX  A 


DERIVATION  OF  ELEMENT  OP  IMPEDANCE  MATRIX 


A.  Derivation  of  Equation  (10) 

From  Equation  (9)  the  impedance  elements  have  the  general  form: 


Z,  =  -  ?  *1^  ds 

In  n 


(A.l) 


where  ?  is  the  expansion  mode  of  the  surface  current  3  and  is  the 
n  s 

radiated  electric  field  generated  by  the  electric  current  mode 
can  be  calculated  by  the  vector  potential  as 


^1  •  VV-A 


=  -jw  (i 


>1 

A  +  - TT- 


Thus,  Equation  (A.l)  becomes 


Zi  =  j 
In 


i»  4’  [' 


^  VV-l^]ds 
Ic 


s 

j»  p^-  %  4  ‘ 

J  k. 


(A.2) 
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By  using  the  vector  identity  one  can  express 

as 

Thus,  Equation  (A. 2)  becomes 

'  p  r 

Zin  =  j«  ^  ^  V*  (?j^V*l^)ds  -  ^  (V*l^)(7*l^)dS. 

,  k  J  J 

(A. 3) 

From  the  surface  form  of  the  Divergence  Theorem,  the  second  term  in 
(A. 3)  reduces  to  the  following  contour  integral: 


^  V- V'^^jdS  =  o  (^^7*l^)*n  dl 

,  k  Jc 


(A. 4) 


where  n  is  the  outward  normal  of  the  integration  contour  C  as  shown  in 
Figure  A.l. 

Because  of  the  symmetry  of  the  geometry,  J^^only  has  t  and  ^ 
components.  For  the  4>  component. 


-n  =  0. 
n 


For  the  t  component. 


^  =  0  at  the  open  aperture  since  jf  is  a  piecewise  sinusoidal 


function.  Consequently, 
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(A.5) 


Thus,  from  Equation  (A. 3), 


1  P 

Z,  =  jw  ^  4^  ds  -  if  (V-1  )(V-^^)ds 
in  n  J  ^ 


(A.6) 


The  vector  potential  A  is  given  by 


OPEN  APERTURE 


Figure  A.l.  Geometry  for  the  surface  form  of  the  Divergence  Theorem. 
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in  which  since  ^j^is  the  current  on  S'  while  7  is  operating  on  S. 

From 
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Equation  (A. 9)  becomes, 


% 

L  1 

> 

V 

L-jkH  1 

7* 

■f  e 

ds'  =  - 

J^-V' 

e 

'^1  R 

^  J 

R 

^  / 

^  -jkR 

S  g . . 

R 


^-jkR 

~R~ 


V'  -J, 


ds' 


(A. 10) 


The  first  term  of  Equation  (A. 10)  can  be  shown  as  zero  by  changing 
the  surface  integral  into  a  contour  integral  and  using  the  same  argument 
as  the  ones  which  lead  to  (A. 5).  Thus, 


-jkR  ) 


1  R 


ds'  = 


,-jkR 

nr 


ds' 


(A.ll) 


and  becomes 


'In  ■  4n 


t  t  n  .n  1 


V'S 


-jkR 

(7'*Jj^)  g-^ —  ds'ds 


4n 


jkR 


^  (V-lj^)(V'-?^)  g-^  ds'ds  (A. 12) 


as  written  in  Equation  (10). 
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B.  Derivation  of  Impedance  Elements 


From  Equation  (9), 


In 


t  ds 

n 


,t<|) 

^In 


<|)  •  E^^  ds 

n 


In 


t  •  E***^  ds 
n 


"In 


<}>  J***  •  e'*’**’  ds 
n 


in  which, 


t  “n<'> 


J*^  =  sin  m<()  P  (t) 


(A. 13) 


and  defined  in  Equations  (6)  and  (7).  By  using  Equation 

(A. 12), 
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,tt 

^In  4ii 


*  I-  *  1  -  t  "  t 


z«+ .  i!!!!”  fi 

^In  4ii  U 


^  (V*<|)jj)(7'-t'jf)  ds'ds 

n  i  ,  Z  n  1  R 


'In  4ll 


*  ^  ds'ds 

Ic 


(A. 14) 


Note  that  ^  ft  ^  P  I? 


Thus, 


’'•‘-'J  -  it)  It  ('‘'’i)  ■  it)  I?  S 


kKt)  "♦  “i<*> 


(A.15) 


’•W  ■  iu  !♦  <■'*>  ■  io  li  "♦  vo) 


p(t)  a* 


=  ^  m<|.  P^(t) 


(A.16) 


Also, 


t  =  x(fp)(p*x)  +  y(t-p)(p-y)  +  z(t*z) 


sin^  cos<|)  +  y  sin^  sin<|)  +  z  cos? 
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where  ^  is  the  angle  between  t  and  z. 


<l>  =  -X  sin<t>  +  y  cos<)). 


Similarly, 


t'  =  X  sin^'  cos*'  +  y  sin^'  sin^'  +  z  cos^' 


<j)  =  -X  sin<j>'  +  y  cos^i' . 


Consequently, 


t*t'  =  sin^  sin^'  cos(i(>-<t>' )  +  cosK  cost,' 


t*<|>'  =  sin^  sin(<{)-()>' ) 


<|)*t'  =  -  sin^'  sin(<(>-(J>' ) 


=  cos(<t>-<ti')  .  (A. 17) 


The  distance  R  between  and  is  given  by 
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¥2 

[2  .  2  21 

(pcos<t>  -  p'cos<|>')  +  (psin<|>  -  p'slni|>')  +  (z-z')  J 

=  j^p^cos^^  +  p'^cos^^'  -  2pp'  cos<}'  cos^' 

+  p^  sin^«|>  +  p'^sin^<|>'  -  2pp'  sin«|)  sin<j)'  +  (z-z')^] 

% 

=  ^p^+p'^  -  2pp'  cos(<|>-4>')  +  (z-z')^j 
=  j^(p^  -  2pp'  +p'^  )  +  (z-z')  +  2 pp' (1-cos («)>-'(•' ))J 
=  [(P-P')^  +  (z-z')^  +  4pp'  sin^  (^^]] 


(A. 18) 


where 


p  =  p(t),  p'  =  p'(t') 
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(1)  Z 


'  ■'^''■'5  -  h  (’■'  ■'91 


=  [t.t'j  cos  ^  cos  mr 


.  2  Uco( 


-jkR 

fcm  "*  ”♦'  “i**'']  TT- 


but  ds'  =  p'(t')  di|)'  dt' 
ds  =  p(t)  d<|>  dt, 


tt  \ 

'in  “  An  lsin^sin^'cos(<(>-<|)' )  +  cos^cos^' jcos  nuj)  cos  in<J)' 

JtJt'J<j>Jr  '■ 


Qjj(t)Q^(t')) 


cos  m<|>  cos  m<|)'  Q^(t)Q'(t')  — —  d(>>'d<}>  dt'dt 


(A. 19) 
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For  the  horn  radiation  problem,  m^.  Consider 

g-jkR 

cos(4>-4*')  cos  m<fr  cosm<|>'  — —  d<|>'  df 

q-3^ 

cos (<!>-♦') [cos  m(<j)-^')  +  cosm(«H-+')]  — ^ d4>' 

(A. 20) 
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(711 


=  n 


e-jkR 
cos\J/cosm\J;  — g —  d\|/ 


2ji 


g-jkR 

cosvj/cosmij/  — —  dij/ 


(|»=0 


=  n 


cos(m+l)»J/  +  cos 


-jkR 


di(/ 


\j;=0 


Let  G  be 
m 


G  = 
m 


e-jkR 

cosin4>  — 5 —  d^* 


thus, 


-jkR 


e  ^ 

cos(4>-<P')  cosm<j>  cosm<()'  — —  d4''d<fr  =  ’^rm+l  ^  ^m-lj 


Consider 


e"^*^  1 

cosm<(>  cosm<(>'  ^  d<}>'d<()  -  2 

r. 

r  1e 

cosm( (()-<{>' )+cosm(^+4>' )  — —  d<f>d<()' 


(A. 22) 


By  letting  !}>=(()-<()',  (A. 22)  becomes 
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1 

2 


[2n  f2ii-^ 


♦'=0j 


cosmif;  + 


cosm(t|H^2<()' ) 


e-jkK 

B 


d\|d^' 


[In 


=  n 


cosm\|> 


d\|/ 


=  2  n  G 

m 


Substituting  Equations  (A. 21)  and  (k.23)  into  (A. 19)  gives 


j«Mo 

4ji 


2cos^cos5' 


Wt'dt 


j«Mo 

■  4k2 


tj 


|^sln^ln^'(G^^j^+G^_j^)  +  2cos^cos^' t)Qj^( 


Wt'dt 


(A. 23) 


Qn(t)Qi(t') 


t') 


(A. 24) 
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(2) 


'‘It 


J‘*'c 

4ji 


n  1 


-jkR 


-72  ('•♦Jn]  ■'1)]  T- 


4ii 


sin  m<f)  Pjj(t)  cos  tn<J)' 


Ql(t') 

kp'(t') 


-jkR 


-sin^' sin(  )  sinm<t>cosm4'' 


P  (t) 
n 


Ql(t') 

kp'(t') 


11  m 

j^2  ^  p(t)p'(t') 


cosra^cos!n(|)' 


Pjj(t)Q£(t')]  • 


-jkR 

p(t)p'(t')d<t.'d<j)dt'dt  (A. 25) 


Consider  (m^) 


sin(4>-4'' )sin  m<}i  cos 


<)>J 


r 


-jkR 


n  r 


1 

2 


jkR 


r  He 

sin( )  sinm(  (()+<f>' )  +  s inm ((})-<)>' )  — = —  d4>d<f> 
L  J  A' 
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Since  ,  the  above  equation  becomes 


sin(<|)-+Osin  m<|>  cos  m(|>' 


<I>J 


-jkR 


[2n 


<|>'=0, 


i  X  2n 


P"-*'  r  . 

sin\|/  sinm(\j;+2(|)' )  +  sinm»|/l  — g —  d«|«i<|)' 


r2ii 


e-jkR 

sini|;  sinmi|;  — —  d\|/ 


<|»=0 


rzn 


i|/=0 


g-jkR 

[-cos(m+l)\|/  +  cos(m-l)»|;]  — ^ —  d\|«J<|)' 


g-jkR 

[cos(m+l)i|»-cos(m-l)\|;]  — g —  dij; 


"[Vl  -  «n,-l] 

(A. 26) 


Substituting  (A. 26)  and  (A. 23)  into  (A. 25)  gives 
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^In  =  An 


tj 


^  G^PJt)Qi(t')]dt'dt 


tJ 


h<')^in5'('=».rViK(')¥'') 


G/„<t)Qi(f)]dfdt 


ik.27) 


(3)  Z 


In' 


■"In  An 


[t  J^-i'  4  -h  jJ)]  dsMs 


3^ 

An 


t'J 


■frj 


r 

|^sin^sin(<j)-i})' )cos  m(}i 


-  ^  bhr  "* 


jkR 


(tTF)  P(t)p'(t')  d4.'d<^  dt'dt 


(A. 28) 


It  can  be  shown  that 
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♦J 


-jkR  /  >. 

sin(<|)-<|>' )  cos  m<|>  sin  m<|>'  —  d<|>'d<|>  =  (A* 29) 


Substituting  (A. 29),  (A. 23)  into  (A. 28)  gives 


^In  An 


tj 


jfE^G„Q'<t)Pj(t')]p(t)p'(t')  dt'dt 


2  kp(t)p'(t') 


J<^c 

4k^ 


tJ 


2m 

k 


G|jjQ;^(t)Pi(t')]dt'dt 


(A. 30) 


(4)  Z 


H. 


In 

Z 


<|)(|>  _  ^***^0 
in  “  4lt 


-jkR 


[♦  ■^1  -  j  ■'i)]  V- 


P  P  P 


4n 


tJ 


<f>J 


|^cos(<|>-<|)' )  sin  m(j»  Pjj(t:)  sin  m<|)'  Pj^(t') 


-  h  yh)  "♦  ‘’n<'>)('7TTO  ’■l*'’’)] 


-jkR 


R 


p(t)p'(t')  d<t>'d<l>  dt'dt  (A. 31) 
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Again,  it  can  be  shown  that 


cos(^-<j>')  sin  m<(>  sin  m4>' 


-jkR  ,  . 


(A.32) 


Substituting  {A.32),  (A. 23)  into  (A. 31)  gives 


Z 


An 


['P(t)‘>'<‘'>(WVl)Pn<t)fl<'') 

tjt' 


2m 


-jkR 


dt'dt 


4k^ 


tj 


[k2p(t)p'(t.)(G,^l*G^_l)-2m\]p„(,)Pi(t')  dt-dt 

(A. 33) 
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